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Abstract 

We present a simple formula for the generating function for the polynomials in 
the d-dimensional semiclassical wave packets. We then use this formula to prove the 
associated Rodrigues formula. 


1 Introduction 

The generating function for 1-diniensional semiclassical wave packets is presented in formula 
(2.47) of [2]. In this paper, we present and prove the d-dimensional analog. As an application 
of this formula, we also prove the associated (multi-dimensional) Rodrigues formula. 

These results have also been proven from a completely different point of view by Helge 
Dietert, Johannes Keller, and Stephanie Troppmann. See Lemma 3 and Section 3 (particu¬ 
larly Proposition 16) and formula (13) of p. See also [3]. We have also received a conjecture 
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from Tomoki Ohsawa |1] that the generating function result could be proved abstractly by 
using the formula for products of Hermite polynomials and the action of the metaplectic 
group. 
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2 Semiclassical Wave Packets 


The semiclassical wave packets depend on two invertible d x d complex matrices A and B 
that are always assumed to satisfy 

A*B + B*A = 21 and A^ B - B^ A = 0. 

They also depend on a phase space point (a, rj) that plays no role in the present work. After 
chosing a branch of the square root, we dehne 


<yCo(A, B, h, a, rj, x) = vr (det A) 

{{x — a), B A~^ {x — a)) . {rj, {x — a)) 


X exp 


2h 


I 


h 


Here, and for the rest of this paper, we regard as being embedded in C'^, and for any two 
vectors a G and b G C'^, we use the notation 


d 

{a, b) = 7Tj bj. 
i=i 


For 1 < / < d, we dehne the raising operator 

TZi = Ai{A, B, h, a, r])* = -^= {{Bci, {x-a)) - i {Aci, {-i KV - r ]))) 
Then recursively, for any multi-index /c, we dehne 


(pk+ei{A, B, h, a, ? 7 , x) = 




TZi^ipkiA, B, h, a, r])){x). 
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For fixed A, B, h, a, rj, these wave packets form an orthonormal basis indexed by k. It 
is easy to see that 

^k{A, B, h, a, r], x) = Pk{A, h, (x - a)) B, h, a, r], x), 

where Pk{A, h, (x — a)) is a polynomial of degree \k\ in (x — a), although from this dehnition, 
it is not immediately obvious that Pk{A, h, (x — a)) is independent of B. 

Since they play no interesting role in what we are doing here, we henceforth assume a = 0 
and rj = 0. 


3 The Generating Function 


Our main result for the generating function is the following: 


Theorem 3.1 The generating function for the family of polynomials Pk{A, h, x) is 
G{x, z) = exp ^ — (z, A~'^ Ax') + (^z, . 


he., 


G{x, z) = Pk{A, h, x) y. 


Remark We make the unconventional dehnition |A| = \JA A* . By our conditions on the 
matrices A and B, this forces |y4| to be real symmetric and strictly positive. We also have 
the polar decomposition A = |A| f/^, where U a is unitary. With this notation, we can write 

2 

\/h 

This equivalent formula is the one we shall actually prove. 

We begin the proof of Theorem 3.1 with a lemma that provides an alternative formula 
for TZi. From this formula and an induction on \k\, one can easily prove that Pk{A, h, x) is 
independent of B, because 

<y9o(A, B, h, 0, 0, x) B, K 0) 0) x) = fi-^P I dety4|“^ exp 


^ {x, |kl| 



(T(x, z) = exp 


Uaz, Uaz) 
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Lemma 3.2 For any ip ^ S 


{Riip){x) 


1 

V 2 B, h, 0, 0, x) 


(^Aei, V ((^o(^, B, h, 0, 0, x) ip{x)) ^ . 


Proof; The gradient on the right hand side of the eqnation in the lemma can act either on 
the '(po or on the ip. So, we get two terms when we compute this: 

\/f BA ~\) + {x, BA ip{x) 

V j=i 

- (Ae,, (V^)(a;))^. 


The second term here is precisely the second term {—i{Aei, i KV) ip{x))), in the 

expression for {Riip){x). So, we need only show the hrst term here equals the hrst term, 

_ {B Cl, x) ip{x), in the expression for {Riip)[x). 

\/2 h 

To do this, we begin by noting that the hrst term here equals 


1 

Iaci, (ej({ej,BA ^ x) + {x, B A ^ej)))ip{x) 




^ CL _ 

IAci, (ej({ej,BA ^x) + {BA ^ Cj, x)) ) ip{x) 


2 ^ ^ 


2\/^ ^ 


a 

(^Aei, (ej(j^ej,BA ^ x) + {B A~^ ej, x)'j^ ip{x) 


.7 = 1 


d 

X] (^Cj (^{ej,BA ^ x) + {ej, {A~^Y B* x)'^^ ip {x) 


-1 


I I , BA ^ + {A-^y B* . , , 

Aei, --- X ) ip{x) 




Because of the relations satished by A and B, BA ^ is (real symmetric) + i (real symmetric). 
So, its conjugate, BA ^ has this same form. Thus, BA ^ equals its transpose, which 
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is {A B*. So, the quantity of interest here equals 


1 

1 

V2h 

1 

1 


{Aeu B* x) 'ip{x) 

{ei, A* B* x) ^(x) 

(ei, B* x) ii{x) 

{Bei, x) i>{x), 


which is what we had to show. H 


Proof of Theorem 3.1: By an induction on |fc|, we prove there is a order polynomial 
Pk in d variables with the following properties: 


• Pk{,A, h, x) = 

• (I?) 


Pk{\A\ ^x/y/h). 

= pk{\A\~^ x/\/h-UAz)G{x, z). 


The result then follows by setting z = 0. We never compute the polynomial pk because it 
may be complicated. We also use the notation pk{A, h, x) = pk{{\A\~^ x/^/h). 


d \ 


k+ei 


We dehne po = 1. Below, we inductively compute G{x, z). For our second 

condition above to hold, Pk+ei must be dehned via the sum of formulas fl3.ip and 03.21) . This 
uniquely defines the polynomial Pk+er 


For k = 0, the result is trivial since Po{A, h, x) = 1. 


For the induction step, it is sufficient to do the following for an arbitrary positive integer 
I < d: 

Assuming we have already proved these for some k, we prove them for the multi-index k + ei. 


To do this, we begin by noting that 
Pk{A, B, h, 0, 0, x) = 


1 

71 


B, K 0, 0))(a;). 
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Also, 


So, 


5, h, 0, 0, x) = 2 1^1/^ (fc!) Pk{A, h, x) ^po{A, B, h, 0, 0, x). 


Tl'"{cpo{A, B, h, 0, 0))(a;) = 2 1^'/^ Pk{A, h, x) ipo{A, B, h, 0, 0, x). 


Thus, when we apply the raising operator, the polynomial Pk{A, h, x) gets changed to 
^k+ei{A, h, x). 

d^G 

Assuming the induction hypothesis, when we differentiate ——^ with respect to zi, the zi 

oz'^ _ 

derivative can act on the G{x, z) or it can act on the pk{\A\~^ x/\/h — Uaz). When it acts 
on the G{x, z), we obtain 


2 ^|/1| ^x/Vh-UA^^ Pk{\A\ xj\Ph-JJAx) G{x, z). 


(3,1) 


From the induction hypotheses, this is a polynomial of degree |/c| + 1 evaluated at the 
argument \A\~^x/-\/h— Uaz. Note that this result depends on the following calculation, 
with G{x, z) written with the polar decomposition of A: 


dzi 


{x, z) 


- {UAeuUAz) - {UAZ,UAei) + -^{UAei, \A\ ^x) ) G{x, z) 


2 (^UAei, (^\A\ ^ x/Vh- UAZ^'j G{x, z). 


When the 
\k\-l. 


acts on the polynomial, we either get zero or a polynomial of degree 

OZi 


- ((Vpfc)(|A| ^ x/y/h-llAz), UaCi'^ G{x, z) 
= - (^UaCi, {Vpk){\A\~^ x/y/h- Uaz)^ G{x, z). 


(3.2) 


Recall that 

{Ri'ip){x) = - -1 /5 ^ (Aei,V (MA B, h, 0, 0, x) V'(x)) ) , 

V 2 B, h, 0, 0, x) \ ^ 2 / 
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and that from our induction hypothesis, 


B, h, 0, 0, x) ^pkiA, B, h, 0, 0, x) 

= 7r-'^/2;j-d/2 2-|fe|/2^^,)-i/2|^g^ A\-^pk{A, h, x) 


exp < — 


{x, \A\ ^ x) 

h 


When computing TZi pk by this formula, the gradient in TZi can act on the exponential or 
the Pk{,A, h, x). When it acts on the exponential, we get 

2 -I'=I /2 (/d)-V2 \A\~'^x) po{A, B, h, 0, 0, x) 

= 2-(l^l+b/2 y^((fc + e^)!)-i/2 

X 2 (^UAeu\A\~^ x/^/Tl^ pk{A, h, x) ipo{A, B, x). (3.3) 

When the gradient in TZi acts on the Pfc(A, h, x), in TZi (pk, we get the term 
- 2"l'=b2 (/j!)-V2 Vx{pk{A, h, x)) ) po{A, B, h, 0, 0, x) 

= - 2 "(l^l+b /2 (/;;!)-V2 /aci, ^ (ej, {Vpk){A, h, x)) \A\-^ eA pq{A, B, h, 0, 0, x) 


i=i 


= - 2 (l^l+b/2 (/^!) 1/2 (^Aeu |A| ^ (Vpfc)(A, h, x) ) ipo{A, B, h, 0, 0, x) 
= - 2-(l'^l+b/2 {{k + 

X {UAei, {\/pk)iA, h, x)) po{A, B, h, 0, 0, x). 


(3.4) 


From fl3.ip and fl3.2p with z = 0, we obtain 


UAei, \A\ ^xjVh) pk{A, h, x) - lUAeu {Vpk){\A\ ^x/Vh) 


From fl3.3p and fl3.4p and taking into account the factor of y/ki + 1 in 
TZi^pk) = \/ki + l(Pk+ei, we obtain 

Pk+ei{A, h, x) 

= 2 (^UAei, \A\~'^ x/Vh^ pk{A, h, x) - l^UAeu iypk){\A\~^x/\fh) 
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The quantities of interest contain the same polynomial evaluated at the appropriate 
arguments, and K = Pk+ei{A, h, x). Since I is arbitrary, with 1 < / < d, the 

result is true for all multi-indices with order \k\ -|- 1, and the induction can proceed. I 


4 The Rodrigues Formula 

As an application of the Generating Function formula, we prove a multi-dimensional Ro¬ 
drigues formula for these polynomials. The result is 

Theorem 4.1 In d-dimensions, with the convention that |A| is the positive square root of 
A A*, 



Remark By scaling, it is sufficient to prove this for h= 1. 

We begin the proof of this result with a lemma that embodies a special case of the chain 
rule for high order derivatives in d-dimensions. 

Lemma 4.2 Assume F : ^ is G™ and that M : 'Mf ^ is linear. Then viewing 

gradients as eolumn vectors and using multi-index notation, 

(VJ” FiMx) = ((M‘ 

One proves this by induction on |m|, proving that each component is correct. 

Proof of Theorem 4.1 Using this technical result and the generating function. 


Pfc(A, 1, x) 

= {V.f e^v(-(U~Azf (TTaz) + 2(t7Jz)‘(|A|-n))^_„ 

= (vyexp(-(f7;z - \A\-'x)‘(jrAZ - |2i|-n)) exp((|2i|-n)‘(|2i|-n))^,„ 
= exp((|/i|-n)‘(|/i|-n)) (vyexp(-(c7ix - |xi|-‘x)‘(rxx - 
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= exp((|A|-‘i)‘(|A|-‘a:)) 

= exp((|xir‘i)‘(|xl|-‘x)) 
= exp((|xir‘i)‘(|xl|-‘x)) 
= exp((|xir‘i)‘(|xl|-‘i)) 

= exp((|xl|-‘i)‘(|xl|-'i)) 


((t'x)‘V„)‘’exp(-(uJ - |xl| 'x)‘{w - |xl| 
((U^A)‘V„)‘’exp(- {w - u)'{w - ti)) 

u = 

(- (?7a)* V„)^exp (-(w - uf{w - u)) ^ = 0 

u = |A|“^x 

{-{U^y\A\^V^)^exp{-{w - \A\-^xy{w - |x4|-A))^^q 

(- (T^y Al* V^)%xp (- (A|-A)* (A|-A)). 


However, \ A\^ = |x4| is real, and with our unusual convention, A = |x4| Ua- So, (I/a)* AI = A*. 
This proves the theorem. I 
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